CTP-TAMU-65/92 
NUB-TH-3055/92 



hep-ph/93023T7| 



COSMOLOGICAL CONSTRAINTS AND SU(5) 
SUPERGRAVITY GRAND UNIFICATION 

^C) ' R. Arnowitt 

On" 

ON ■ Center for Theoretical Physics, Department of Physics 

^ : Texas A&M University, College Station, TX 77843-4242 



IX, 
(N 



X 



and 



Pran Nath 



ABSTRACT 



' Department of Physics, Northeastern University 

Boston, MA 02115 

m 

(N 
O 

m 

The predictions of SU(5) supergravity models with radiative breaking constrained 

Oh! 

by experimental proton decay bounds are discussed. It is shown that cosmological 
constraints further restrict the parameter space but can be satisfied for a wide range 
of parameters. It is also shown that no serious fine tuning problems (either at Msusy 
or Mqut) exist. 



1. INTRODUCTION 

The observation last year [1] that the SU(3)c x SU(2)l x U(1)y gauge coupling con- 
stants, Q!3, «2 and a i = (5/3)ay, meet at a common energy scale fx = Mq if extended 
from their measured values at \i = Mz by the supersymmetric renormalization group 
equations (RGE) with one pair of Higgs doublets, has led to a number of investigations 
of other predictions of supersymmetric GUT models [2-9]. In supergravity grand unifica- 
tion [10] with radiative breaking of SU(2) x U(l) [11] (for reviews see [12]), the model 
depends on four parameters (aside from the as yet unknown t-quark mass mj): uiq (uni- 
versal scalar mass), mi/2 (universal gaugino mass), Aq (cubic soft breaking mass) and 
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tan (3 =< H 2 > / < H 1 >, where < H 2 > gives mass to the up quarks and < H\ > to 
the down quarks and leptons. Thus with gauge group SU(5), Refs. [5 and 3] discuss the 
No-scale model Aq = = mo, Ref. [7] the case Bo = Ao — mo (where So is the quadratic 
soft breaking mass and can be expressed in terms of tan f3) and Refs. [2,4,6,7] examine 
the general parameter space. Ref. [9] is concerned with the O (10) model. 

The fact that supergravity grand unification introduces only four unknown parameters 
(only two more than in the Standard Model itself) to account for the masses of 32 particles 
(31 new SUSY particles plus the light Higgs h) implies that there should be considerable 
correlation between the SUSY masses. Unfortunately if one imposes only the requirement 
that radiative breaking occur and the current experimental bounds on the SUSY masses, 
the allowed mass bands are still very broad [2,5,7] and it is difficult to make clear predictions 
that can be used to test the theory. The situation changes considerably, however, if the 
model possesses an <S77(5)-type proton decay: p — > vK. Thus if one assumes no extreme 
fine tuning of parameters (mo, rrig < 1 TeV where g is the gluino) and the superheavy 
Higgs color triplet, which mediates the decay obeys M# 3 < 3Mq (which in simple models 
is the bound that keeps the GUT couplings perturbative in size), then the parameter 
space allowed by current data is still fairly large e.g.: mo ^ 550 GeV, m~ g ^ 450 GeV 
(i.e. mi/2 ^ 150 GeV), 1.1 < tan/3 < 4.7 [2]. However, one finds a number of remarkable 
predictions for the SUSY masses [2]: 

2m 2l ^ m 22 m^ ~ (- - -)m § 
m z 3 ~ m z 4 ~ m w 2 

where the charginos (Wj, i = 1,2) and neutralinos (Zi, i = 1 . . .4) are labeled such that 
nii < rrij for i < j. In addition one finds rrih ^ 110 GeV and m t ^ 180 GeV. Further, for 
m t < 140 GeV, then m^ ^ 100 GeV when m^ ^ 95 GeV making one of these particles 
(and possibly both) observable at LEP200. 

The above proton decay constraint is sufficiently powerful to eliminate the preferred 
models of Ref. [7] and the No-scale SU(5) model over the entire parameter space [3]. 
(For the No-scale case one need not even use the fine tuning constraint if one imposes the 
cosmological requirement that the LSP be electrically neutral.) Thus No-scale models are 
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viable only if they can suppress the p — > vK decay mode, such as is done in the flipped 
SU(5) supergravity model [13]. 

The purpose of this letter is to discuss the role of cosmology in a GUT theory which 
allows proton decay via dimension five operators. We shall show that a supergravity 
GUT theory which is constrained both by proton decay limit and the cosmological relic 
density limit (which avoids overdosing the universe) allows a wide domain of the parameter 
space on a reduced four dimensional manifold (more precisely a five dimensional shell) in 
contradiction to the conclusions of a recent analysis on this topic [14]. We also show that 
the conclusions drawn in Ref. [14] concerning fine tuning are inaccurate, and there are no 
serious fine tuning problems either at Msusy or Mq. 

2. COSMOLOGICAL CONSTRAINTS 

Recently detailed analyses of the neutralino relic density in N=l supergravity unified 
models have been carried out using the superparticle spectrum generated by the radia- 
tive electroweak symmetry breaking [14-16]. It is found that the relic density constraint 
^Zi^ 2 ~ 1 ( w here Oj, is the ratio of the lightest neutralino mass density to the critical 
mass density and h is Hubble constant measured in units of 100 km/sec Mpc) limits ad- 
ditionally the allowed parameter space of the supergravity models when one imposes the 
naturalness condition discussed above, m^, mo ^ 1 TeV. [If squark and gluino masses in 
excess of 5 TeV are allowed, the relic density constraint is easily satisfied [15] (due to the 
many open channels), as is the p-decay constraint (due to the suppression from the large 
SUSY masses in the dressing loop [17]).] 

A detailed analysis of the allowed parameter space under the simultaneous constraints 
of proton stability [17] and neutralino relic density not overdosing the universe [18] will be 
given elsewhere [19]. We give here a brief discussion. For an arbitrary point in parameter 
space one finds O^r « 100. As has been pointed out [14-16], however, the neutralino 
annihilation rate is significantly enhanced when the annihilation occurs close (within a few 
GeV) to the h boson being on shell, i.e. 2m^ ~ m^. To calculate the Z\ relic density 
near the h pole, we follow the general analysis of Ref. [18] making use of the cross section 
for Zi + Zi — > h* —>/ + / where / is a final state fermion [20]. However, near the 
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pole, it is necessary to take the thermal average of the rigorous cross section, < av > 
(irrelative velocity), as discussed in detail in Ref. [21], rather than use the expansion 
av = a + bv 2 /6. The thermal average can no longer be done analytically, but must be 
performed numerically. (To our knowledge, previous calculations of SUSY relic densities 
have not included this important modification.) From Eq. (1), the condition that the 
intermediate h is nearly on-shell can be viewed as a constraint relating m g to m^, and we 
find strong suppression of fi^/i 2 over a range of gluino masses ^ 5-20 GeV wide. Thus 
the inclusion of the relic density constraint reduces the five dimensional parameter space, 
mo, nig, Aq, tan [3 and rrit to a four dimensional submanifold (actually a five dimensional 
shell ?s 5-20 GeV wide) where the annihilation of the Z\ is enhanced so that fi^/i 2 < 1- 
If m t is experimentally determined (as one hopes it soon will be at the Tevatron) then one 
will be left with a three dimensional subspace (more precisely a four dimensional shell) 
depending on the parameters mo, Aq, and tan (3. 

While the relic density constraint reduces the allowed range of the remaining parame- 
ters somewhat, significantly it still leaves available a wide range of these parameters. Fig. 
1 shows the allowed region in the parameters m g , an (tan (3 = l/tana#) for a charac- 
teristic example mo = 600 GeV, A t = 0, fx > where A t is the t-quark A parameter at 
the electroweak scale. (In Figs. 1 and 2, we have allowed the more conservative bound of 
M H3 < 6 Mq.) We see that m g ranges from 200 GeV to 450 GeV and 22° < a H < 41°, 
much as when only the p-decay constraint was imposed [2]. (The allowed m g band is ?s 10 
GeV wide.) Note also that allowed ranges of parameters satisfying both relic density and 
proton decay constraints exist for a full range of values of m t and not for only "special" 
values of m t as stated in Ref. [14] . Fig. 2 shows the dependence of the allowed region 
on m g and A t for mo = 600 GeV and an = 30° (tan (3 = 1.73). The m g band is again 
10-20 GeV wide and the range of A t is similar to that obtained before [2]. Again, allowed 
regions in parameter space exist for a range of values of m t . Fig. [3] shows the importance 
of using the rigorous analysis of Ref. [21] rather than the approximation < av >= a + bx 
(x = kT '/m^ )■ As one can see from comparison of the exact and the approximate analysis, 
the approximate analysis would introduce large errors. 

The above analysis shows that the relic density constraint combined with the proton 
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decay constraint for arbitrary m t leads to a wide range of allowed parameters, in contrast to 
the analysis in Ref. [14]. Further, the model has many experimentally testable predictions. 
There is, however, an alternate framework which retains all the predictions of the standard 
SU(5) model in terrestrial experiments, and eliminates the constraint of the neutralino relic 
density altogether by a sufficiently rapid decay of the neutralinos. We begin by recalling 
that we are dealing with an N = 1 supergravity theory which is an effective remnant theory 
at the scale Mq of some more unified structure. As such, it can possess at the scale Mq 
operators with dim > 3 in the superpotential scaled by Mpg or the compactification mass 
Mc- While the effect of these operators on the computation of the superparticle spectrum 
via radiative breaking [2-9] would be negligible, they can significantly affect the cosmology 
resulting from the model. 

As an example, consider the second generation of quarks and leptons and supplement 
the particle content of the model by an SU(5), singlet field v c and an SO (10) singlet 
field N so that they can be grouped into nonets of SU(3)c x SU(3)l x SU(3)r. We 
denote these nonets by the representation L(l, 3, 3), Q(3, 3, 1) and Q c (3, 1, 3) of [SU(3)] 3 . 
Next we extend the superpotential of the theory by adding the minimal terms W s + 
A[(Tr£ 2 ) / 'Mo\TrQLQ c , where W s contains the singlet fields v c and N and generates super- 
heavy VEVs for them, and E is the 24 of 5*7(5) whose VEV, diag< E >= M (2, 2, 2, -3, -3) 
breaks SU(5) to SU(3)c x SU(2)l x U(1)y at Mq- The second term contains the factor 

TrQLQ c = —DND C + D£ c u c - Du c d c - qlD c + qHu c + qH'd c (2) 

where D, D c are the superheavy Higgs color triplets, I and q are lepton and quark fields, 
and H and H' are the two light Higgs doublets. We note that both SU(5) and [.St/^)] 3 
can be embedded into Eq, pointing to a possible common origin of the normal and "Planck 
slop" terms from a more unified structure. 

After spontaneous breaking, W s generates a VEV for v c which spontaneously produces 
a violation of R-parity as can be seen from Eq. (2). This also gives superheavy masses 
to the N and v c fields. The v c VEV growth generates a D — d and D c — d c mixing. 
Diagonalization of the D — d, etc. mass terms [22] gives an effective R-parity violating 
interaction in the quark- lepton sector which can decay the lightest neutralino. Thus after 
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diagonalization, the R-parity violating interaction determining this decay is 



C int = g(Z^ L )(s a cl) ■ g = — f 22 U Sl (3) 

m a 

where fi(x) is the muon field, s a (x), c a (x) (a = color index) the strange and charm quark 
fields, f 22 the SU{5) coupling in f l3 H x M lY Mf Y (H x = 5 Higgs field, M Y = 5, M XY = 10 
matter fields, i,j = generation indices), U is the projection of the Z\ state onto the photino 
state, si is the D — d mixing parameter given by 

s 1 = -y/(m 2 D + y 2 ) 1 / 2 ; V = ~^2~ < v c > \ (4) 

and mo = M Ha is the superheavy Higgs color triplet mass. The partial lifetime of the Z\ 
is then 



r(Z x - asp*) ~ (1 x 10- 19 sec)(/ 22Sl t/)- 2 (^) 4 (^^) (5) 

with f 22 = (m s /M z )(e/ sma H sm28 w ). For < E >=< v c >= M G ^ 10 16 GeV, M c = 
5 x 10 17 GeV, M D = 3M G , A = 1, U = 1, = 500 GeV, m 2l = 50 GeV, one finds 
t « 10~ 4 sec. Thus typically the neutralino is very short lived so that it will not leave any 
significant cosmological trace. However, the lifetime of the neutralino is still large enough 
that it will decay well outside the detection chamber in collider experiments, so that all of 
the characteristic missing Et signals of supersymmetric particles will be maintained. 

3. FINE TUNING AT M SU sy 

The problem of fine tuning first arose in non-SUSY GUTs due to the quadratic di- 
vergence of ma, the Higgs mass i.e. m 2 H = m\ — bah 2 , where mo is the bare mass, 
a is a coupling constant, A is the cutoff and b is a constant. Thus if mjj = O (Mz), 
then one must fine tune m\ to 24 decimal places when A = Mq. One may formal- 
ize this argument [23,14] by defining the parameter c = {m^/m 2 H ){dm 2 H /dmo). Then 
c = m\jm 2 H = cak 2 /m 2 H ~ 10 24 , i.e. log c is the number of fine tuning decimal places. In 
supersymmetry the fine tuning problem resurfaces in the radiative breaking equation [11] 
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X -Ml = (m 2 Hl - m 2 H2 tan 2 /?)/(tan 2 /3 - 1) - . (6) 

Here ji is the H\ — H 2 mixing mass, mH 1>2 are the Higgs masses which can be expressed 
in terms of the GUT scale parameters mo, m^, A®, Bq, /iq by the renormalization group 
equations. The situation here is more complex as there are many parameters = m 2 ,, /Uq, 
my 2 etc - One ma y define Cj = (ai/M^)(dM^/dai) and require Cj < Aj with, say Aj = 10 2 
as in Ref. [14]. There are however, a number of ambiguities that need to be addressed. 
Thus one can always make transformations on the parameters, a\ = fi(cii) sending a to 
c[, increasing or decreasing the value of a given a in this way. (The in general are 
complicated functions of the hidden sector of the theory [12]. Thus which set of functions 
of <2j are "fundamental" and hence to be preferred is unknown at present.) Also rescaling 
all parameters to a single one e.g. by writing £o = nio/mi/2, £a = ^-0/^1/2 etc. (as 

done in Ref. [14]) artificially increases the remaining c-parameter as Ci/ 2 then equals E q. 

In Ref. [14], the choice = \i 2 and m\ was made. The authors then found only c t 
"too high" (and then only by a factor of 2-3). We believe it is incorrect to use m t as a 
fine tuning parameter as (presumably) the top will shortly be discovered, and then the 
only correct thing would be to insert its experimental value. But even allowing the choice 
a = m 2 , one could reduce the value of c t by a factor of two merely by replacing M\ by 
Mz in the definition of c t , which would then satisfy the criteria of Ref. [14]. 

The above discussion shows that the fine tuning criteria used in Ref. [14] is ambiguous 
up to factors of 2-10. When one fine tunes to 24 decimal places as in non-SUSY GUTs, these 
ambiguities are unimportant. But if one is talking about conditions such as A < 10 2 , no 
clear conclusions can be drawn. The only reasonable constraint is the qualitative one used 
in Ref. [2] that squark and gluino masses be less than 1 TeV (which is also approximately 
the detection upper bound at the LHC and SSC). 

4. FINE TUNING AT M G 

We mention briefly some additional points concerning the problem of fine tuning at 
the GUT scale. There are three theoretically satisfactory methods of breaking SU(5) to 
the Standard Model group while maintaining light Higgs doublets and superheavy Higgs 
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triplets. The first is that originally proposed in global SUSY GUT models [24], which 
requires a fine tuning that, however, is natural due to the SUSY no renormalization the- 
orems. These models use a VEV growth of a 24 representation to break SU(5). The 
second is the missing partner models [25] using 50, 50 and 75 representations of SU(5) 
where the VEV growth in the 75 breaks 577(5). The third method makes use of a global 
577(6) symmetry in the GUT sector [26]. The breaking of the GUT group then makes the 
Higgs doublets pseudo Goldstone bosons and hence automatically light. We view this last 
method as being more elegant than the missing partner models in either flipped SU(5) or 
normal 577(5). 

The third method illustrates the fact that fine tuning in the GUT sector may not 
be as invidious as other fine tunings. Thus the physics of the GUT sector is at present 
unknown, and some higher symmetry (perhaps from string theory) may naturally force 
two coupling constants to be equal, thus keeping the Higgs doublets light: this is precisely 
what happens in the SU (6) model above. 

5. DISCUSSION 

At present there is no acceptable string model, Calabi-Yau, orbifold or 4-D construc- 
tion, which is consistent with the coupling constant unification analysis [1]. Thus the 
fact that a GUT model may possibly have string anticedents appears irrelevant at this 
point, as it may have to be significantly modified if and when a viable string model is 
constructed. More relevant is the "possibility" that the No-scale model could determine 
the soft breaking parameters dynamically. For the complete No-scale model, where only 
the mi/2 soft breaking mass is non-zero at the GUT scale, this would lead to a unique 
SUSY mass spectrum at a fixed value of m t . Implementing this is the most interesting 
question facing the No-scale model, since then very likely it could be experimentally deter- 
mined whether it is right or wrong. Barring this theoretical development, one should look 
for experimental differences between different supergravity models. One of these involves 
proton decay, where the flipped SU(5) model suppresses the p — > VK decay [13], while 
SU (5) supergravity expects it to be seen at the Super Kamiokande experiment. 

This work was supported in part by NSF Grant Nos. PHY-916593 and PHY-917809. 

8 



Note added: After completing this work we received the preprint Ref. [27]. These authors 
also find a significant region of parameter space simultaneously satisfying the SU(5) proton 
decay and cosmo logical constraints. They also reconfirm the mass relations of Eq. (1). 

FIGURE CAPTIONS 

Fig. 1. Region in m g — an parameter space allowed by the combined proton decay and Z\ 
relic density constraints for mo = 600 GeV, A t = 0.0, > (tan (3 = I/tana^)- M is 
the t-quark A parameter at the electroweak scale. The dashed curve is for m t = 110 
GeV, solid curve for m t = 125 GeV, dot-dash for m t = 140 GeV. 

Fig. 2. Region in m g — A t parameter space allowed by the combined proton decay and Z\ 
relic density constraints for mo = 600 GeV, an = 30° (tan (3 = 1.73), fx > 0. Different 
curves as in Fig. 1. 

Fig. 3. Contribution to fi^/i 2 from the h pole for m t = 125 GeV, m = 600 GeV, A t /m = 
0.5, tan (3 = 1.73, /j > 0. The solid line is the exact result from thermal averaging 
over the Higgs pole. The dashed line is the approximate result when one expands 
av ^ a + bv 2 /6. 
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